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In an analysis of a 2.92 fb 1 data sample taken at 3.773 GeV with the BESIII detector operated at 
the BEPCII collider, we measure the absolute decay branching fractions to be B(D° —>■ K~e + v„) — 
(3.505 ± 0.014 ± 0.033)% and B(D° n~e + u e ) = (0.295 ± 0.004 ± 0.003)%. From a study of the 
differential decay rates we obtain the products of hadronic form factor and the magnitude of the 
CKM matrix element f+{0)\V cs \ = 0.7172±0.0025±0.0035 and fl(0)\V cd \ = 0.1435±0.0018±0.0009. 
Combining these products with the values of |V^ s (rf) | from the SM constraint fit, we extract the 


3 


hadronic form factors f+{ 0) = 0.7368 ± 0.0026 ± 0.0036 and /+(0) = 0.6372 ± 0.0080 ± 0.0044, 
and their ratio /+(0)// + (0) = 0.8649 ± 0.0112 ± 0.0073. These form factors and their ratio are 
used to test unquenched Lattice QCD calculations of the form factors and a light cone sum rule 
(LCSR) calculation of their ratio. The measured value of /^^(OjlWsCd)! and the lattice QCD value 
for /|' (7r) (0) are used to extract values of the CKM matrix elements of || = 0.9601 ± 0.0033 ± 
0.0047 ± 0.0239 and \V c d\ = 0.2155 ± 0.0027 ± 0.0014 ± 0.0094, where the third errors are due to the 
uncertainties in lattice QCD calculations of the form factors. Using the LCSR value for f+(0)/f+ (0), 
we determine the ratio | |/1 V^ s | = 0.238 ± 0.004 ± 0.002 ± 0.011, where the third error is from 

the uncertainty in the LCSR normalization. In addition, we measure form factor parameters for 
three different theoretical models that describe the weak hadronic charged currents for these two 
semileptonic decays. All of these measurements are the most precise to date. 

PACS numbers: 13.20.Fc, 12.15.Hh 


I. INTRODUCTION 

In the Standard Model (SM) of particle physics, the 
mixing between the quark flavors in the weak interaction 
is parameterized by the unitary 3x3 Cabibbo-Kobayashi- 
Maskawa (CKM) matrix Vckm [1, 2]. The CKM matrix 
elements are fundamental parameters of the SM, which 
have to be measured in experiments. Beyond the SM, 
some New Physics (NP) effects would also be involved 
in the weak interactions of the quark flavors, and mod¬ 
ify the coupling strength of the quark flavor transitions. 
Due to these two reasons, precise measurements of the 
CKM matrix elements are very important for many tests 
of the SM and searches for NP beyond the SM. Each 
CKM matrix element can be extracted from measure¬ 
ments of different processes supplemented by theoretical 
calculations of corresponding hadronic matrix elements. 
Since the effects of the strong and weak interactions can 
be well separated in semileptonic D° —> K~e + u e and 
D° —y ir~e + v e decays, these processes are well suited for 
the determination of the magnitudes of the CKM matrix 
elements V cs and V c d, and also for studies of the weak 
decay mechanisms of charmed mesons. If any significant 
inconsistency between the precise direct measurements 
of \V c d\ or |V CS | and those obtained from the SM global 
fit is observed, it may indicate that some NP effects are 
involved in the first two quark generations [3]. 

In the limit of zero positron mass, the differential rate 
for D° —> K~(-K~)e + v e decay is given by 

^ wi»-(.->i 3 i/f ( A di 2 , im 

where Gp is the Fermi coupling constant, Pk-(-k~) is 
the three-momentum of the meson in the rest 

frame of the D° meson, and (q 2 ) represents the 

hadronic form factors of the hadronic weak current that 
depend on the square of the four-momentum transfer 
q = pp,o — Pk~(tt-)- These form factors describe strong 
interaction effects that can be calculated in lattice quan¬ 
tum chromodynamics (LQCD). 

In recent years, LQCD has provided calculations of 
these form factors with steadily increasing precision. 
With these improvements in precision, experimental val¬ 


idation of the computed results are more and more im¬ 
portant. At present, the main uncertainty of the apex 
of the Bd unitarity triangle (UT) of B meson decays is 
dominated by the theoretical errors in the LQCD de¬ 
terminations of the B meson decay constants /b ( s) and 
decay form factors /+^ 7r (0) [3]. Precision measurements 

of the charmed-sector form factors fp^\q 2 ) can be used 
to establish the level of reliability of LQCD calculations 
of /f^(0). If the LQCD calculations of fp^\q 2 ) agree 

well with measured f+^\q 2 ) values, the LQCD calcula¬ 
tions of the form factors for B meson semileptonic decays 
can be more confidently used to improve measurements 
of B meson semileptonic decay rates. The improved mea¬ 
surements of B meson semileptonic decay rates would, in 
turn, improve the determination of the Bd unitarity tri¬ 
angle, with which one can more precisely test the SM and 
search for NP. 

In this paper, we present direct measurements of the 
absolute branching fractions for D° —> K~e + v e and 
D° —► 7r _ e + ^ e decays using a 2.92 fb -1 data sample 
taken at 3.773 GeV with the BESIII detector [4] oper¬ 
ated at the upgraded Beijing Electron Positron Collider 
(BEPCII) [5] during the time period from 2010 to 2011. 
(Throughout this paper, the inclusion of charge conjugate 
channels is implied.) By analyzing partial decay rates for 
D° —> K~e + v e and D° —> n~e + v e , we obtain the q 2 de¬ 
pendence of the form factors /^^(q 2 ). Furthermore we 
extract the form factors (0) and /+ (0) using values of 
|T4s| and |T4d| determined by the CKMfitter group [6]. 
Conversely, taking LQCD values for f+(0) and /+(0) as 
inputs, we determine the values of the CKM matrix ele¬ 
ments \V CS \ and | V cd \■ 

We review the approaches for describing the dynamics 
of D° —> K~e + v e and D° —> Ti~e + v e decays in Section II. 
We then describe the BESIII detector, the data sample 
and the simulated Monte Carlo events used in this analy¬ 
sis in Section III. In Section IV, we introduce the analysis 
technique used to identify the semileptonic decay events. 
The measurements of the absolute branching fractions for 
these two decays and study of systematic uncertainties in 
these branching fraction measurements are described in 
Section V. In Section VI, we describe the analysis tech¬ 
niques for measuring the differential decay rates for these 
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two semileptonic decays, and present our measurements 
of the hadronic form factors. The determinations of the 
CKM matrix elements \V CS \ and | V c ,i | are discussed in 
Section VII. We give a summary of our measurements in 
Section VIII. 


II. FORM FACTOR AND APPROACHES FOR 
D° SEMILEPTONIC DECAYS 


(KS) [ 8 ] and Bauer-Stech-Wirbel (BSW) [9] models, a 
commonly used form factor has a single pole of the form 

f + (q 2 ) = /+( °J , (II.3) 

1 

which is simply the first term in Eq. (II.2) (taking a = 0). 
The pole mass M po \ e is often treated as a free parameter 
to improve fit quality. 


A. Hadronic form factor 


In general, the form factor f+^\q 2 ) can be expressed 
in terms of a dispersion relation [7] 


/+(<7 2 ) = 


/+(0)/(l-a) | 1 

1 - 


^7 


h+ t-q 2 -ie 


(II 1) 


where M po \ e is the mass of the lowest-lying relevant vec¬ 
tor meson, for D° —> K~e + v e it is the D* + , while for 
D° —> n~e + p e it is the D* + , /p(0) is the form factor 
evaluated at the four-momentum transfer q = 0 , a is 
the relative size of the contribution to /+( 0 ) from the 
vector pole at q 2 = 0 , t+ = (mjjo + nriK-(-R-)) corre¬ 
sponds to the threshold for D°K~(n~) production, rrijj o 
and ’rrix-(n-) are the masses of the D° and charged kaon 
(pion) meson, respectively. From Eq. (II. 1) we find that, 
except for the pole position of the lowest-lying meson be¬ 
ing located below threshold, f+(q 2 ) is analytic outside of 
a cut in the complex < 7 2 -plane extending along the real 
axis from t + to oo, corresponding to the production re¬ 
gion for the states with the appropriate quantum num¬ 
bers. 


B. Parameterizations of form factor 

The form of the dispersion relation given in Eq. (II.1) 
is often parameterized by keeping the lowest-lying me¬ 
son pole explicitly and approximating the remaining dis¬ 
persion integral in Eq. (II. 1) by a number of effective 
poles [7, 10] 


/+( 9 2 ) = 


/ + (0)/(l-q) , 1 
1 - 




N 

■£ 

k =i 


Pk 


i - 


(II.2) 


where p^ and 7 *, are expansion parameters that are not 
predicted. The form factor can be approximated by intro¬ 
ducing arbitrarily many effective poles. Equation (II.2) 
is the starting point for many proposed form factor pa¬ 
rameterizations. 


1. Single pole form 


2. Modified pole model 


The modified pole model uses N = 1 in Eq. (H-2), and 
the form factor can be expressed as 


/+(<? 2 ) 


/+( 0 ) 




Q 


(II.4) 


where a is a free parameter. This model is the so-called 
Becirevic-Kaidalov (BK) parameterization [10] and has 
been used in many recent lattice calculations and exper¬ 
imental studies for D° meson semileptonic decays. 


3. Series expansion 


The series expansion [7] is the most general parame¬ 
terization that is consistent with constraints from QCD. 
It has the form 


/+(*) = to ) ao{to) + > 

(II.5) 

where 


Z{t, T 0 = —==- . 

\/t+ ~ t + V4 — 1 

t- = ( m D o - m K -( n -)) 

t 0 = 4(1 — a /1 — *-/*- 


(11.6) 

(11.7) 

(11.8) 


a 0 (t 0 ) and Tk{t 0 ) are real coefficients. The function 
Pit) = z{t,m 2 D ») for D K and P{t) = 1 for D —> tt. 
$ is given by 


$(Mo) = 


1 




-5 


24t txv t+ ~ to 


X (\A+ ~t+ y/t+ — to){y/t+ — t + \/t+ — t _) 3 ^ 2 

x (4 - tf/\ (II.9) 


where \v can be obtained from dispersion relations us¬ 
ing perturbative QCD and depends on the ratio of the 
s quark mass to the c quark mass, £ = m s /m c [11]. At 
leading order, with £ = 0 , 


In the constituent quark model, lattice gauge calcula¬ 
tions, and QCD sum rules, such as the Korner-Schuler 


3 

327r 2 m 2 


( 11 . 10 ) 
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The choice of P and <f> is such that 

a o(^o) ( 1 + Y2 r fe(^) ) < !■ (11.11) 


fc=1 


The 2 series expansion is model independent and satisfies 
analyticity and unitarity. In heavy quark effective the¬ 
ory [12] the coefficients in Eq. (II. 5) for D —> ire + v e 
and B —> ir£ + ve decays are related. A measurement of 
the Tk for the decay of D —> tt e + v e therefore provides im¬ 
portant information to constrain the class of form factors 
needed to fit the decays of B —> 7r£ + ^, and thereby pro¬ 
vides improvements in the determination of the magni¬ 
tude of the CKM matrix element V u b . However, the valid¬ 
ity of the form factor parameterization given in Eq. (II. 5) 
still needs to be checked with experimental data. This is 
one of the reasons why it is important to precisely mea¬ 
sure the form factors /+ ^(q 2 ) for D° —> K~{-K~)e + v e 
decays. 

In practical applications, one often takes fc max — 1 or 
fcmax = 2 in Eq. (II. 5), which gives following two forms 
of the form factor: 

(a) Series expansion with 2 parameters of the form 
factor is given by 


/+W P(t)$(Mo) a ° (<o) 

x (l + ri(io)[z(t,i 0 )]), (11.12) 


which gives 


f (t) = 1 / + (0)P(0)$((Mo) 

+ P(t)$(t,t Q ) 1 + ri(t 0 )z(0,i 0 ) 

x (l + ri(t 0 )[z(Mo)]). (11.13) 


(b) Series expansion with 3 parameters of the form 
factor is given by 


/+(i) = P(t)$(Mo) a ° (io) 

x (l + n(to)[z(t,t 0 )] 

+ r 2 (t 0 )[z(t,t 0 )] 2 ), (H.14) 


which gives 


P(t)$(t,t 0 ) 

x / + (0)P(0)4>(0,to) 

1 + ri(t o )z{0,to) +r 2 (to)[-(0,to)] 2 
x (1-Hr!(t 0 Mi,to)] 

+r 2 {t 0 )[z(t,t 0 )] 2 ). (11.15) 


III. DATA SAMPLE AND THE BESIII 
EXPERIMENT 

At y/s = 3.773 GeV, the i/>(3770) resonance is di¬ 
rectly produced via e + e _ annihilation. About 93% [6] of 
i/j( 3770) decays to DD ( D°D °, D + D~) meson pairs. In 
addition, the continuum processes e + e _ —t qq (q = u, d , s 
quark), e + e _ — > t + t~ , e + e _ — > 7isr J/4>, e + e~ — > 
7isr' ! /’(3686) events are also produced, where 7isr is the 
radiative photon in the initial state. The data sample 
contains a mixture of all these classes of events. In the 
analysis, we refer to events other than ^>(3770) decays to 
DD as “non -DD process” events. 

BEPCII [5] is a double-ring e + e~ collider operat¬ 
ing in the center-of-mass energy region between 2.0 
and 4.6 GeV. Its design luminosity at 3.78 GeV is 
10 33 cm _2 s _1 with a beam current of 0.93 A. The peak 
luminosity of the machine reached 0.65 x 10 33 cm” 2 s _1 at 
y/s = 3.773 GeV in April 2011 during the ^(3770) data 
taking. BESIII [4] is a general purpose detector operated 
at the BEPCII. At the BEPCII colliding point, the e + 
and e~ beams collide with a crossing angle of 22 mrad. 

The BESIII detector is a cylindrical magnetic detec¬ 
tor with a solid angle coverage of 93% of 47r. It con¬ 
sists of several main components. Surrounding the beam 
pipe, there is a 43-layer main drift chamber (MDC) that 
provides precise measurements of charged particle tra¬ 
jectories and ionization energy losses ( dE/dx ) that are 
used for particle identification. The momentum resolu¬ 
tion for charged particles at 1 GeV/c is 0.5%, and the 
specific dE/dx resolution is 6%. Outside of the MDC, 
a time-of-flight (TOF) system is used for charged par¬ 
ticle identification. The TOF consists of a barrel part 
made of two layers with 88 pieces of 2.4 m long plas¬ 
tic scintillators in each layer, and two end-caps with 96 
fan-shaped detectors. The TOF time resolution is 80 ps 
in the barrel, and 110 ps in the end-caps, corresponding 
to a K/ 7r separation better than 2cr for momenta up to 

1 GeV/c. An electromagnetic calorimeter (EMC) sur¬ 
rounds the TOF and is made of 6240 CsI(Tl) crystals ar¬ 
ranged in a cylindrical shape (barrel) plus two end-caps. 
The EMC is used to measure the energies of photons and 
electrons. For 1.0 GeV photons, the energy resolution is 
2.5% in the barrel and 5.0% in the end-caps, and the one¬ 
dimensional position resolution is 6 mm in the barrel and 
9 mm in the end-caps. A superconducting solenoid mag¬ 
net outside the EMC provides a 1 T magnetic field in the 
central tracking region of the detector. A muon identifi¬ 
cation system is placed outside of the detector, consisting 
of about 1272 m 2 of resistive plate chambers arranged in 
9 layers in the barrel and 8 layers in the end-caps incor¬ 
porated in the magnetic flux return iron of the magnetic. 
The position resolution of the muon chambers is about 

2 cm. This system efficiently identifies muons with mo¬ 
mentum greater than 500 MeV/c over 88% of the total 
solid angle. 

The BESIII detector response was studied using Monte 
Carlo event samples generated with a GEANT4-based [17] 
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detector simulation software package, BOOST [18]. To 
match the data, 1.98 x 10 8 Monte Carlo events for 
e + e _ — > t/j( 3770) — > DD were simulated with the Monte 
Carlo event generator /C/C, KKMC [19], where 56% of the 
ip( 3770) resonance is set to decay to D°D° while the 
remainder decays to D + D~ meson pairs. All of these 
D°D° and D + D~ meson pairs are set to decay into differ¬ 
ent final states which were generated with EvtGen [20] 
with branching fractions from the Particle Data Group 
(PDG) [ 6 ]. This Monte Carlo event sample corresponds 
to about 11 times the luminosity of real data. With these 
Monte Carlo events, we determine the event selection cri¬ 
teria for the data analysis and study possible background 
events for the measurement of the D° —»■ K~e + v e and 
D° —> T:~e + v e decays. We refer to these Monte Carlo 
events as “cocktail vs. cocktail DD process” events. 

Since the “non -DD process” events are mixed with 
the DD events in the data sample, we also generate 
“non -DD process” Monte Carlo events simulated with 
KKMC [19] and EvtGen [20] to estimate the number of 
the background events in the selected D° — > K~e + v e and 
D° —> ir~e + v e samples. 


A. Properties of doubly tagged D°D° decays 

For a specific tag decay mode, the number of the single 
D° tags is given by 

Atag = 27V£,0£)o/?tagetagj (IV.1) 

where N d oqo is the number of the D°D° meson pairs 
produced in the data sample, Stag is the branching frac¬ 
tion for the tag mode, and e tag is the efficiency for recon¬ 
struction of this mode. Similarly, the number of the D° 
semileptonic decay events observed in the system recoil¬ 
ing against the single D° tags is given by 

A observed —> h e + V e ) = 2N D oQo BtagB(D° — > h e + u e ) 

^ Uag,D° — >h~ e + is e i (IA.2) 

where h~ denotes the final state hadron (i.e. h~ = K~ 
or 7T“), B(D° —> h~e + v e ) is the D° meson semileptonic 
decay branching fraction, and e tag ,D°^h~e+u e is the ef¬ 
ficiency of simultaneously reconstructing both the single 
D° tag and the D° meson semileptonic decay. With these 
two equations, we obtain 


To estimate the efficiencies, we also generate “Signal” 
Monte Carlo events, i.e. ip( 3770) —> D°D° events in 
which the D° meson decays to all possible final states [ 6 ] , 
and the D° meson decays to a semileptonic or a hadronic 
decay final state that is being investigated. These Monte 
Carlo events were all generated and simulated with the 
software packages mentioned above. 


IV. RECONSTRUCTION OF D°D° DECAY 
EVENTS 


B(D° -*• h~e + v e ) 


A 0 bserved ( D® > h e~t^ e ) , o\ 

N tag e(D°^h-e+v e ) ’ 1 j 


where e(D° ->■ h e + v e ) = e tag , D o^ h -e+vj £ta g - 

To measure the D° semileptonic differential decay rates 
given in Eq.(I.l) we need to evaluate the partial decay 
rate AT,; observed within a small range of the squared 
four-momentum transfer A qf, where i stands for the ith 
q 2 bin. This partial decay rate can be evaluated with the 
double tag D°D° events as well. The measurement of 
the partial decay rates is described in Section VI. 


B. Single D° tags and efficiencies 


In ^>(3770) resonance decays into DD mesons in which 
a D meson is fully reconstructed, all of the remaining 
tracks and photons in the event must originate from the 
accompanying D. In these cases, the reconstructed me¬ 
son is called a single D tag. Using the single D° tag 
sample, the decays of D° —> K~e + v e and D° —> n~e + p e 
can be reliably identified from the recoiling tracks in the 
event. We refer to the event in which the 5° meson is re¬ 
constructed and a semileptonic D° decay is reconstructed 
from the recoiling tracks as a doubly tagged D°D° decay 
event or a double D°D° tag. With these doubly tagged 
D°D° events, the absolute branching fractions and the 
differential decay rates for D° semileptonic decays can 
be well measured. 

In the analysis, all 4-momentum vectors measured in 
the laboratory frame are boosted to the e + e _ center-of- 
mass frame. 

In this section, we describe the procedure for selecting 
the single D° tags and the D° semileptonic decay events. 


The D° meson is reconstructed in five hadronic decay 
modes: K + tt~, K + n~Tr°, , A' + 7 r _ 7 r + 7 r _ 7 r° 

and Ar + 7 r _ 7 r° 7 r°. Events that contain at least two re¬ 
constructed charged tracks with good helix fits are se¬ 
lected. The charged tracks used in the single tag analy¬ 
sis are required to satisfy |cos0| < 0.93, where 6 is the 
polar angle of the charged track. All of these charged 
tracks are required to originate from the interaction re¬ 
gion with a distance of closest approach in the trans¬ 
verse plane that is less than 1.0 cm and less than 15.0 cm 
along the z axis. The dE/dx and TOF measurements 
are combined to form confidence levels for pion (CL 
and kaon ( CLk ) particle identification hypotheses. In 
the selection of single D° tags, pion (kaon) identifica¬ 
tion requires CL V > CLk ( CLk > CL for momenta 
p < 0.75 GeV/c and CL V > 0.1% ( CL K > 0.1%) for 
p > 0.75 GeV/c. 

A 7T° meson is reconstructed via the decay 7r° —> 77 . 
To select photons from n° decays, we require an energy 
deposit in the barrel (end-cap) EMC to be greater than 
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0.025 (0.050) GeV and in-time coincidence with the beam 
crossing. In addition, the angle between the photon and 
the nearest charged track is required to be greater than 
10°. A one-constraint (1-C) kinematic fit is performed 
to constrain the invariant mass of 77 to the mass of 7 r° 
meson, and \ 2 < 50 is required. 

For the D° —> K + tt~ final state, we reduce back¬ 
grounds from cosmic rays, Bhabha and dimuon events 
by requiring the difference of the time-of-flight of the two 
charged tracks be less than 5 ns, and the opening angle 
of the two charged track directions be less than 176 de¬ 
gree. In addition, we require that the sum of the ratio of 
energy over the momentum of the charged track is less 
than 1.4. 

The single D° tags are selected using beam energy con¬ 
strained mass of the Kmr (where n = 1,2, 3, or 4) com¬ 
bination, which is given by 

Mbc = (IV.4) 

where Flbeam is the beam energy and \pKmr\ is the mag¬ 
nitude of the momentum of the daughter Kmr system. 
We also use the variable A E = Exun — -Fbeam, where 
Eru 7 t is the energy of the Kmr combination computed 
with the identified charged species. Each Kmr combina¬ 
tion is subjected to a requirement of energy conservation 
with |A.E| < (2 ~ 3)&e k „„, where crE Kn „ is the standard 
deviation of the Ektut distribution. For each event, there 
may be several different charged track (or both charged 
track and neutral cluster) combinations for each of the 
five single D° tag modes. If more than one combination 
satisfies the energy requirement, the combination with 
the smallest value of |AE| is retained. 

The dots with error bars in Fig. 1 show the resulting 
distributions of Mbc for the five single D° tag modes, 
where the D° meson signals are evident. To determine 
the number of the single D° tags that are reconstructed 
for each mode, we fit a signal function plus a back¬ 
ground shape to these distributions. For the fit, we use 
signal shapes obtained from simulation convolved with 
a double-Gaussian function for the signal component, 
added to an ARGUS function multiplied by a third-order 
polynomial function [ 21 , 22 ] to represent the combinato¬ 
rial background shape. The ARGUS function is [23] 


/argus ( m ) 


= m 




1 - 



(IV. 5) 

where m is the beam energy constrained mass, Eq is the 
endpoint given by the beam energy and c is a free param¬ 
eter. The solid lines in Fig. 1 show the best fits, while 
the dashed lines show the fitted background shapes. 

In addition to the combinatorial background, there 
are also small wrong-sign (WS) peaking backgrounds 
in single D° tags. The doubly Cabibbo suppressed 
decays (DCSD) contribute to the WS peaking back¬ 
ground for single D° tag modes of D° —> K + n ~, D° —> 
K + Tr~TT° and D° —> K + ir~Tr~Tr + . In addition, the 


D° KgK~ir + (K% 7 T+ 7 T-), D° K° S K~ tt+tt 0 

(Kg —* 7r + 7r - ) and KgK~ir + (Kg —>• 7r°7r°) also make 
significant contributions to WS peaking backgrounds for 
the D° —>■ K + ir~ tt~ tt + , D° K + ir~Tr~Tr + Tr° and 
D° —> K + Tr~TT°Tr° tag modes, respectively. The size of 
these WS peaking backgrounds are estimated from Monte 
Carlo simulation and then subtracted from the yields ob¬ 
tained from the fits to Mbc spectra. 

Table I summarizes the single D° tags. In the table, 
the second column gives the A E requirement on the Kmr 
combination, the fourth column gives the number of the 
single D° tags in the tag mass region as shown in the 
third column. 

The efficiencies for reconstruction of the single D° 
tags for the five tag modes are obtained by applying 
the identical analysis procedure to simulated “Signal” 
Monte Carlo events mixed with “Background” Monte 
Carlo events. The “Signal” Monte Carlo events are gener¬ 
ated as e + e _ —» '0(3770) —¥ D°D° , where the D° meson 
is set to decay to the tag mode in question and the D° 
meson is set to decay to all possible final states with cor¬ 
responding branching fractions [ 6 ]. The efficiencies for 
reconstruction of the single D° tags are presented in the 
last column of Table I. 



M bc (Gev/c 2 ) 


FIG. 1. Distributions of the beam energy constrained masses 
of the Kmr (n = 1, 2, 3 or 4) combinations for the 5 single 
D° tag modes: (a) K + ir~, (b) A + 7r _ 7r°, (c) A + 7r _ 7r _ 7r + , (d) 
A + 7r _ 7r _ 7r + 7r° and (e) A' + 7r _ 7r°7r°. 


C. Selection of D° —y K e^~h> e and D° —> ir e+iz e 

The D° —> K~e + v e and D° —> Tr~e + v e event can¬ 
didates are selected from the tracks recoiling against 





















TABLE I. Summary of the single D° tags and efficiencies for reconstruction of the single D° tags, where A E gives the 
requirements on the energy difference between the measured Ekuix and beam energy Ebeam, while the Mb range defines the 
signal region of the single D° tags. A r tag is the number of single D° tags and etag is the efficiency for reconstruction of the 
single D° tags. 


Tag mode 

A E (GeV) 

Mbc range (GeV/c 2 ) 

Wag 

«tag (%) 

K+ 7T~ 

(-0.049, 0.044) 

(1.860,1.875) 

567083 ± 848 

70.29 ± 0.07 

K + TT~Tr° 

(-0.071,0.052) 

(1.858,1.875) 

1094081 ± 1692 

36.80 ± 0.03 

K + 7T _ 7T _ 7T + 

(-0.043,0.043) 

(1.860,1.875) 

700061 ± 1121 

39.57 ± 0.04 

A' + 7T _ 7r _ 7r + 7r° 

(-0.067,0.066) 

(1.858,1.875) 

158367 ± 749 

15.95 ±0.08 

O 

o 

1 

+ 

(-0.082,0.050) 

(1.858,1.875) 

273725 ± 2859 

15.78 ±0.08 

Sum 



2793317 ± 3684 



the single D° tags. To select the D° —> K e + u e and 
D° —>• 7r e + v e events, it is required that there are only 
two oppositely charged tracks, one of which is identified 
as a positron and the other as a kaon or a pion. The 
combined confidence level CLk ( CL n ) for the K (ir) 
hypothesis is required to be greater than CL n r (CLk) 
for kaon (pion) candidates. For positron identification, 
the combined confidence level ( CL e ) : calculated for the 
e hypothesis using the dE/dx , TOF and EMC measure¬ 
ments (deposited energy and shape of the electromag¬ 
netic shower), is required to be greater than 0.1%, and 
the ratio CL e /(CL e ± CL V ± CL K ) is required to be 
greater than 0.8. We include the 4-momenta of near-by 
photons with the direction of the positron momentum to 
partially account for final-state-radiation energy losses 
(FSR recovery). In addition, to suppress fake photon 
background it is required that the maximum energy of 
any unused photon in the recoil system, E 7jmax , be less 
than 300 MeV. 

Since the neutrino escapes detection, the kinematic 
variable 


f^miss — ^miss Prnissj (IV.6) 

is used to obtain the information about the missing neu¬ 
trino, where E m j ss and p m i ss are, respectively, the total 
missing energy and momentum in the event, computed 
from 


Euiiss — F'beaTfi Ijh - Ef ; + , (IV. 7) 

where E h - and E e + are the measured energies of the 
hadron and the positron, respectively. The p m i ss is cal¬ 
culated by 


Pmiss = Pd° - Ph- ~ Pe+ , (IV.8) 

where Pd°, Ph- an< I Pe+ are the momenta of the D° me¬ 
son, the hadron and the positron, respectively. The 3- 
momentum p^o of the D° meson is computed by 

Pd° = -Pta g ^/-E^ eam - m 2 D0 , (IV.9) 

where p tag is the direction of the momentum of the single 
D° tag. If the daughter particles from a semileptonic 



U miss (GeV) 


FIG. 2. U miss distributions of events for (a) D° tags vs. 
D° —> K~e + i'e, and for (b) D° tags vs. D° —>■ Tr~e + u e , 
where the dots with error bars show the data, the solid lines 
show the best fit to the data, and the dashed lines show the 
background shapes estimated by analyzing the “cocktail vs. 
cocktail DD process” Monte Carlo events and the “non -DD 
process” Monte Carlo events (see text for more details). 


decay are correctly identified, C/miss is zero, since only 
one neutrino is missing. 

Figures 2 (a) and (b) show the U m i ss distributions for 
the D° —> K~e + v e and D° —> -n~e + v e candidate events, 
respectively. In both cases, most of the events are from 
the D a —y K~e + v e and D° —> Tr~e + i , e decays. Back¬ 
grounds from DD processes include mistagged D° and 
D° decays other than the semileptonic decay in ques¬ 
tion. Other backgrounds are from “non-Hi) process” 
processes. From the simulated “cocktail vs. cocktail DD 
process” events, we find that the DD background events 
are mostly from D° —> K~n 0 e + i' e , D° —> K~p 1 + u fl 
and D° —> n~e + iy e selected as D° —> K~e + is e , and 
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D° —> n~TT 0 e + v e , D° -A K~e + v^ and D° —> 7 t~ fi + v^ 
selected as D° —> n~e + v e . Backgrounds from “non -DD" 
processes include the ISR (Initial State Radiation) return 
to the 0(3686) and J/0, continuum light hadron produc¬ 
tion, -0(3770) —► non — DD decays and e + e“ -A t + t~ 
events. The levels of these backgrounds events are es¬ 
timated by analyzing the corresponding simulated event 
samples. 

Because of ISR and FSR (Final State Radiation), the 
signal U m i ss distributions are not Gaussian; instead, the 
I/miss distributions have Gaussian cores with long tails at 
both the lower and the higher sides of the distributions. 
To obtain the numbers of the signal events for these two 
semileptonic decays, we fit these distributions with an 
empirical function that includes these tails. 

We use the same probability density function as 
CLEO [24] for t/ miss , 

! ai (of ^ a i + 2 ) if x > ai 

exp(— x 2 /2) if — a 2 < x < cti, 

a 2 (Sf - “2 - a) if x < -a 2 

(IV.10) 

where x = (U m i ss — m)/cr, m and a are the mean value 
and standard deviation of the Gaussian distribution, re¬ 
spectively. In Eq. (IV.10), ai = (m/ai) ni e _ “i/ 2 , and 
a 2 = (n 2 /a 2 ) n2 e _ “ 2 / 2 , where ou, a 2i n\ and n 2 are pa¬ 
rameters describing the tails of the signal function, de¬ 
termined from fits to the simulated t/ m ; ss distributions of 
signal Monte Carlo events. 

To account for differences between data and Monte 
Carlo, we fit the data using the Monte Carlo determined 
f(x) distribution convolved with a Gaussian function 
with free mean and width. The background function is 
formed from histograms of C7 m i ss distributions for back¬ 
ground events from the “cocktail vs. cocktail” DD and 
“non-DZ)” simulated event samples. The normalizations 
of the signal and background are free parameters in the 
fits to the data. 

The results of the fits to the two U m i ss distributions are 
shown in Figs. 2 (a) and (b); the fitted yields of signal 
events are 

^observed (D° -A K~e + v e ) = 70727.0 ± 278.3 (IV.ll) 
and 

^observed(.D 0 -A TT~e + v e ) = 6297.1 ± 86.8. (IV.12) 

In Fig. 2 (a) and (b), the solid lines show the best fits to 
the data, while the dashed lines show the background. 

To gain confidence in the quality of the Monte Carlo 
simulation, we examine the momentum distributions of 
the kaon, the pion and the positron as well as cosdwe 
from the semileptonic decays of D° —> K~e + v e and 
D° —> TT~e + u e , where Owe is the angle between the di¬ 
rection of the virtual W + boson in the D° rest frame 
and the three-momentum of the positron in the W + rest 
frame. These distributions are shown in Figs. 3 (a)-(f), 


respectively, where the dots with error bars are for the 
data, the solid histograms are for the full Monte Carlo 
simulation and the shaded histograms show the Monte 
Carlo simulated backgrounds only. 


V. MEASUREMENTS OF ABSOLUTE DECAY 
BRANCHING FRACTIONS 

A. Efficiency for reconstruction of semileptonic 
decays 

To determine the efficiency e(D° -A h~e + u e ) for re¬ 
construction of each of the two semileptonic decays for 
each single tag mode, “Signal” Monte Carlo event sam¬ 
ples of 0(3770) -A D°D° decays, where the D° meson 
is set to decay to the h~e + i , e final state in question and 
the D° meson is set to decay to each of the five single 5° 
tag modes, are generated and simulated with the BESIII 
software package. By subjecting these simulated events 
to the same requirements that are applied to the data we 
obtain the reconstruction efficiencies £tas,D°^h-e+u e f° r 
simultaneously finding the D° meson semileptonic decay 
and the single D° tag in the same event; these are given 
in Tab. II. 

Due to their low multiplicity, it is usually easier to 
reconstruct D° tags in semileptonic events than in typical 
D°D° events (tag bias). In addition, the size of the tag 
bias is correlated with the multiplicity of the tag mode. 
In consequence the overall efficiencies shown in Tab. II 
vary greatly from the D° -a K~tt + mode to the D° —> 
K~tt + tt + tt~ and D° -a K~n + TT°n 0 modes. 

The last row in Tab. II gives the overall efficiency which 
is obtained by weighting the individual efficiencies for 
each of the five single D° tags by the corresponding yield 
shown in Tab. I. 

There are small differences in efficiencies for find¬ 
ing a charged particle and for identifying the type of 
the charged particle between the data and Monte Carlo 
events that are discussed below in Section VC. To take 
these differences into account, the overall efficiencies 
£mc(D° -a K~e + v e ) and cmc (D° -a TT~e + i / e ) are cor¬ 
rected by the multiplicative factors of 

f trk+pro / 1-0118 for D° -a K~e + v e , , , 

■' corr ( 0.9814 for D° -A 7r _ e + ^ e - ' ’ 

After making these corrections, we obtain the “true” 
overall efficiencies for reconstruction of these two semilep¬ 
tonic decays, 

e(D° -A K~e + v e ) = 0.7224 ± 0.0012, (V.2) 

and 

(V.3) 


e(D° -A Tt-e+ve) = 0.7643 ± 0.0013. 
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FIG. 3. Distributions of particle momenta and cos Owe from D° —¥ K~e + u e and D° —> n~e + i / e semileptonic decays, where 
(a) and (b) are the momenta of kaon and positron from D° —► K~e + i/ e , respectively, (d) and (e) are the momenta of pion and 
positron from D° — > n~e + i/ e , respectively; (c) and (f) are the distributions of cosdwe for D° —> K~e + v e and D° — > n~e + i / e , 
respectively; these events satisfy —0.06 < Umiss < 0.06 GeV. The solid histograms are Monte Carlo simulated signal plus 
background; the shaded histograms are Monte Carlo simulated background only. 


TABLE II. Double tag efficiencies for reconstruction of “Dt ag vs. D° —» h e + v s " and overall efficiencies for reconstruction of 
D° —> h~e + u e in the recoil side of D° tags. 


Tag mode 

^tage+ u e 

^tag,D°— >tt~ e+ v e 

£mc(D° K~e+u e ) 

£mc (D° -¥ n~e + i/ e ) 

K + n~ 

K + n~n° 

K + n~ n~ 

K + n~ n~ n° 
K + 'k~tv°tv 0 

0.4566 ± 0.0014 

0.2685 ± 0.0006 

0.2666 ± 0.0008 

0.1260 ±0.0008 

0.1331 ±0.0007 

0.4995 ± 0.0014 
0.2927 ±0.0007 
0.2897 ±0.0008 
0.1363 ± 0.0008 
0.1467 ±0.0007 

0.6496 ± 0.0021 
0.7296 ± 0.0017 
0.6737 ± 0.0021 
0.7900 ± 0.0064 
0.8435 ± 0.0062 

0.7106 ±0.0021 
0.7954 ± 0.0020 
0.7321 ± 0.0022 
0.8545 ± 0.0066 
0.9297 ± 0.0065 

Average 



0.7140 ±0.0012 

0.7788 ± 0.0013 

B. 

Decay branching fraction 

C. 

Systematic uncertainties 

in measured branching 


fractions 

Inserting the number of the single D° tags, the num¬ 
bers of the signal events for these two D° semileptonic Table III lists the sources of the systematic uncertain- 
decays observed in the recoil of the single D° tags to- ties in the measured semileptonic branching fractions, 

gether with corresponding efficiency into Eq.(IV.3), we We discuss each of these sources in the following, 

obtain the absolute decay branching fractions 


B{D° K~e + v e ) = (3.505 ± 0.014 ± 0.033)% (V.4) 

and 

B(D° t T~e + is e ) = (0.295 ± 0.004 ± 0.003)%, (V.5) 

where the first errors are statistical and the second sys¬ 
tematic. The sources of systematic uncertainties in the 
measured decay branching fractions are discussed in the 
next subsection. 


1. Uncertainty in number of D° tags 

To estimate the uncertainty in the number of single 
D° tags, we repeat the fits to the Mbc distributions by 
varying the bin size, fit range and background functions. 
We also investigate the contribution arising from possible 
differences in the 7r° fake rates between data and Monte 
Carlo simulation. Finally, we assign a systematic uncer¬ 
tainty of 0.5% to the number of D° tags. 
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2. Uncertainty in tracking efficiency 

The uncertainties for finding a charged track are es¬ 
timated by comparing the efficiencies for reconstructing 
the positron, kaon and pion in data and Monte Carlo 
events. 

Using radiative Bhabha scattering events selected from 
the data and simulated radiative Bhabha scattering 
events, we measure the difference in efficiencies for find¬ 
ing a positron between data and simulation. Considering 
both the cos0, where 9 is the polar angle of the positron, 
and momentum distributions of the positrons, we obtain 
two-dimensional weighted-average efficiency differences 
(edata/eMC - 1) of (0.22 ± 0.19)% and (0.11 ± 0.15)%. 
These translate uncertainties on the decay branching 
fractions of 0.19% and 0.15% for D° —> K~e + v e and 
D° —> Tr~e + v e decays, respectively. 

The efficiencies for finding a charged kaon and a 
charged pion are determined by analyzing doubly tagged 
DD decay events. In the selection of the doubly 
tagged DD decay events, we exclude one charged kaon 
or one charged pion track and examine the variable 
Mj( |iss K or n , defined as the difference between the miss¬ 
ing energy squared U r 2 niss and the missing momentum 
squared pjb gs of the selected DD decay events. By ana¬ 
lyzing these MA SS k or T variables for both the data and 
the simulated “cocktail vs. cocktail DD process” Monte 
Carlo events, we find the differences in efficiencies for re¬ 
constructing a charged kaon or a charged pion between 
the data and the Monte Carlo events as a function of the 
charged particle momentum. Considering the momen¬ 
tum distributions of the kaon and pion from these two 
semileptonic decays, we obtain the magnitudes of sys¬ 
tematic differences and their uncertainties of the track 
reconstruction efficiencies. The level of uncertainties in 
the corrections for these differences in measurements of 


TABLE III. Sources of the systematic uncertainties in the 
measured branching fractions for D° —» K~e + v e and D° —> 
n~e + u e . 


Systematic uncertainty (%) 

Source K~e + v e n~e + v £ 


Number of D u tags 

0.50 

0.50 

Tracking for e + 

0.19 

0.15 

Tracking for K~ 

0.42 

— 

Tracking for 7r~ 

— 

0.28 

PID for e+ 

0.16 

0.14 

PID for K~ 

0.10 

— 

PID for 7r" 

— 

0.19 

£" 7 ,max Cut 

0.10 

0.10 

Fit to Umiss 

0.48 

0.50 

Form factor structure 

0.10 

0.10 

FSR recovery 

0.30 

0.30 

Finite MC statistics 

0.17 

0.17 

Single tag cancelation 

0.12 

0.12 

Total 

0.94 

0.90 


the decay branching fractions and partial decay rates (see 
Section VI) are 0.42% and 0.28% for charged kaons and 
pions, respectively. 

3. Uncertainty in particle identification 

The differences in efficiencies for identifying a positron 
between the data and the Monte Carlo samples depend 
not only on the momentum of the positron, but also on 
cos 9. Considering both of these for our signal positrons, 
we obtain a weighted-average difference in efficiency for 
identifying the positron from the two semileptonic de¬ 
cays. After making correction for these differences in effi¬ 
ciencies for identifying the positrons, we obtain a system¬ 
atic uncertainty of 0.16% (0.14%) on the K~(TT~)e + v e 
mode from this source. 

The systematic uncertainties associated with the ef¬ 
ficiencies for identifying a charged kaon and a charged 
pion are estimated using the missing mass square tech¬ 
niques discussed above. Taking into account the mo¬ 
mentum distributions of the charged particles from the 
two semileptonic modes, we correct for the momentum- 
weighted efficiency differences for identifying the kaon 
and the pion, and we assign systematic uncertainties of 
0.10% and 0.19% for charged kaons and pions, respec¬ 
tively. 

4- Uncertainty in E 7jmax cut 

The uncertainty associated with the E imax require¬ 
ment on the events is estimated by analyzing doubly 
tagged DD events with hadronic decay modes. With 
these events, we examine the fake photons from the EMC 
measurements. By analyzing these selected samples from 
both the data and the simulated Monte Carlo events, 
we find that the magnitude of difference in the number 
of fake photons between the data and the Monte Carlo 
events is 0.10%, which is set as the systematic uncertainty 
due to this source. 


5. Uncertainty in fit to U m i SS distribution 

To estimate the systematic uncertainty in the numbers 
of signal events due to the fit to the f7 m i ss distribution, 
we vary the bin size and the tail parameters of the signal 
function. We then repeat the fits to the U m i ss distribu¬ 
tions, and combine the changes in the yields in quadra¬ 
ture to obtain the systematic uncertainty. Since the back¬ 
ground function is formed from many background modes 
with fixed relative normalizations, we also vary the rel¬ 
ative contributions of several of the largest background 
modes based on the uncertainties in their branching frac¬ 
tions and the uncertainties in the rates of misidentifying 
a hadron (muon) as an electron. Finally we find that the 
relative sizes of this systematic uncertainty are 0.48% 








12 


and 0.50% for D° —>K e + v e and D° —> n e + u ei re¬ 
spectively. 

6 . Uncertainty in form factors 

In order to estimate the systematic uncertainty associ¬ 
ated with the form factor used to generate signal events 
in the Monte Carlo simulation, we re-weight the signal 
Monte Carlo events so that their q 2 distributions match 
the measured spectra. We then re-measure the branching 
fraction (partial decay rates in different q 2 bins) with the 
new weighted efficiency (efficiency matrix). The maxi¬ 
mum relative changes in branching fraction (partial decay 
rates in different q 2 bins) is 0.05%. To be conservative, 
we assign a relative systematic uncertainty of 0.10% to 
the branching fraction measurements for D° —> K ~e+v e 
and D° —> n~e + i , e decays. 

7. Uncertainty in FSR recovery 

The difference between the measured branching frac¬ 
tion obtained with the FSR recovery of the positron mo¬ 
mentum and the one obtained without the FSR recovery 
is assigned as the most conservative systematic uncer¬ 
tainty due to FSR recovery. We find the magnitude of 
this difference to be 0.30% for both D° —> K~e + v e and 
D° —> 7T e + v e decays. 

8 . Uncertainty due to finite Monte Carlo statistics 

The uncertainties associated with the finite Monte 
Carlo statistics are 0.17% for both D° —> K~e + v e and 
D° —> n~e + i , e . 


9. Uncertainty due to single tag cancelation 

Most of the systematic uncertainties arising from 
the selection of single D° tags are canceled due to 
the double tag technique. The un-canceled systematic 
error of MDC tracking, particle identification and 
7 r° selection in single tag selection is estimated by 

(Stag( £ tag/ e tag ~ 1) x 0.25<5tag X -^tag) / (^tag ^t a s) ’ 
where e( ag and et a g are the efficiencies of reconstructing 
single D° tags obtained by analyzing the Monte Carlo 
events of D° —>■ tag vs. D° —y h~e + v e and D° —> tag 
vs. D° —> anything after mixing all the simulated 
backgrounds, respectively; N tag is the number of single 
D° tags reconstructed in data; d tag is the total system¬ 
atic error of MDC tracking, particle identification and 
7T° selection in single tag selection. Since no efficiency 
correction is made in the single tag selection, the 
uncertainty in MDC tracking (or particle identification) 
for charged kaon or pion is taken to be 1.0% per track, 


and the uncertainty in 7r° selection is taken to be 2.0% 
per 7T°. For each single 5° tag mode, the uncertainty 
in MDC tracking, particle identification or 7r° selection 
are added linearly separately, and then they are added 
in quadrature to obtain the total systematic error in 
the single D° tag selection. Finally, we assign a sys¬ 
tematic uncertainty of 0.12% for the branching fraction 
measurements. 


D. Comparison with other measurements 

A comparison of our measured branching fractions for 
D° —> K~e + v e and D° —> n~e + i , e decays with those 
previously measured by the MARK-III [25], CLEO [26], 
BES-II [21], CLEO [24, 27, 28] (at the CLEO-c experi¬ 
ment) and BABAR [30, 31] Collaborations as well as the 
world average given by the PDG [6] is given in Ta¬ 
ble IV. Our measured branching fractions for these two 
decays are in excellent agreement with the experimen¬ 
tal results obtained by other experiments, but are more 
precise. In the table, we also compare our branching 
fraction measurements to theoretical predictions for these 
two semileptonic decays. The precision of our measured 
branching fractions are much higher than those of the 
LQCD [14, 32], the QCD sum rule [33] and the LCSR [34] 
predictions. 


VI. DIFFERENTIAL DECAY RATES 

The differential decay rate dr / dq 2 for D° —> 
K~(n~)e + n e is given by Eq. (1.1). The form fac¬ 
tor f+^\q 2 ) can be extracted from measurements of 
dr /dq 2 . Such measurements are obtained from the event 
rates in bins of q 2 ranging from q?—0.5Aq 2 to q 2 +0.5Aq 2 , 
where A q 2 is the bin width and i is the bin number. 


A. Measurement of differential decay rates 

The q 2 value is given by 

q 2 = ( E e + E v ) 2 - (Pe+P,) 2 , (VI.l) 

where E e and p e are the measured energy and momentum 
of the positron, E v and p v are the energy and momentum 
of the missing neutrino: 

Ev = Emiss, (VI. 2) 

Pv = AjnissPniiss- (VI. 3) 

For the D° —> K~e + i' e differential rate, we divide the 
candidates for the decays into 18 q 2 bins. For the D° —> 
n~e + v e mode, which has fewer events, we use 14 q 2 bins. 
The first columns of Tables V and VI give the range 
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TABLE IV. Comparison of the measured B(D° —¥ K e + v e ) and B(D° —¥ n e + v e ) values with those measured by other 
experiments and theoretical predictions based on QCD and the world-average value of r D o. 


Experiment/Theory 

B{D° -A K~e + v e ) (%) 

B{D° ->■ t T~e+u e ) (%) 

PDG2014 [6] 

3.55 ±0.05 

0.289 ± 0.008 

MARK-III [25] 

3.4 ±0.5 ±0.4 

0.391^? ± 0.04 

CLEO [26] 

3.82 ±0.11 ±0.25 

BES-II [21] 

3.82 ± 0.40 ± 0.27 

0.33 ±0.13 ±0.03 

CLEO-c [28] 

3.50 ± 0.03 ± 0.04 

0.288 ± 0.008 ± 0.003 

Belle [29] 

3.45 ±0.07 ±0.20 

0.255 ±0.019 ±0.016 

BABAR [30, 31] 

3.522 ± 0.027 ± 0.045 ± 0.065 

0.2770 ± 0.0068 ± 0.0092 ± 0.0037 

BESIII (this experiment) 

3.505 ±0.014 ±0.033 

0.295 ± 0.004 ± 0.003 

LQCD [14] 

3.77 ±0.29 ±0.74 

0.316 ± 0.025 ± 0.070 

LQCD [32] 

2.99 ± 0.45 ± 0.74 

0.24 ±0.06 

QCD SR [33] 

2.7 ±0.6 


LCSR [34] 

3.9 ± 1.2 

0.30 ± 0.09 





Umiss (GeV) 


FIG. 4. Distributions of I/miss for D° tags vs. D° —> K~e + u e with the squared 4-momentum transfer q 2 filled in different q 2 
bins. The dots with error bars show the data, the blue solid lines show the best fits to the data, while the red dashed lines 
show the background shapes. 
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FIG. 5. Distributions of [/miss for _D° tags vs. D° tv~e + v e with the squared 4-momentum transfer q 2 filled in different q 2 
bins. The dots with error bars show the data, the blue solid lines show the best fits to the data, while the red dashed lines 
show the background shapes. 


of each q 2 bin for D° —> K e + v e and D° —> tt e + v e , 


TABLE V. Summary of the range of each q 2 bin, the num¬ 
ber of the observed events A^bserved, the number of produced 
events N plo duced, and the partial decay rate AT in each q 2 bin 
for D° —> K~e + v e decays. 

q 2 (GeV 2 /c 4 ) A/pbserved_ N produced AF (ns -1 ) 


(0.0, 

0.1) 

7876.1 

± 

94.2 

10094.9 

± 

132.3 

8.812 

± 

0.116 

(0.1, 

0.2) 

7504.3 

± 

90.5 

10015.4 

± 

140.8 

8.743 

± 

0.123 

(0.2, 

0.3) 

6940.5 

± 

87.2 

9502.6 

± 

142.0 

8.295 

± 

0.124 

(0.3, 

0.4) 

6376.0 

± 

83.4 

8667.9 

± 

138.6 

7.567 

± 

0.121 

(0.4, 

0.5) 

6139.8 

± 

81.9 

8575.9 

± 

137.7 

7.486 

± 

0.120 

(0.5, 

0.6) 

5460.5 

± 

77.1 

7384.0 

± 

128.1 

6.446 

± 

0.112 

(0.6, 

0.7) 

5120.3 

± 

74.7 

7101.8 

± 

125.8 

6.200 

± 

0.110 

(0.7, 

0.8) 

4545.5 

± 

70.5 

6322.2 

± 

120.2 

5.519 

± 

0.105 

(0.8, 

0.9) 

4159.4 

± 

67.1 

5760.3 

± 

113.3 

5.028 

± 

0.099 

(0.9, 

1.0) 

3680.7 

± 

63.2 

5183.5 

± 

107.6 

4.525 

± 

0.094 

(1.0, 

1.1) 

3199.6 

± 

58.9 

4550.0 

± 

100.2 

3.972 

± 

0.087 

(1.1, 

1.2) 

2637.1 

± 

53.5 

3810.2 

± 

92.4 

3.326 

± 

0.081 

(1.2, 

1.3) 

2239.1 

± 

49.4 

3239.1 

± 

84.3 

2.828 

± 

0.074 

(1.3, 

1.4) 

1752.1 

± 

43.9 

2621.2 

± 

77.3 

2.288 

± 

0.067 

(1.4, 

1.5) 

1301.0 

± 

37.7 

1989.4 

± 

67.4 

1.737 

± 

0.059 

(1.5, 

1.6) 

927.5 ±32.0 

1505.1 

± 

59.0 

1.314 

± 

0.052 

(1.6, 

1.7) 

541.3 ± 24.6 

983.4 

± 

50.3 

0.858 

± 

0.044 

(1.7, 

(/max/ 

188.2 ± 15.1 

434.2 

± 

39.6 

0.379 

± 

0.035 


respectively. 

The points with error bars in Figs. 4 and 5 show the 
t^miss distributions for the D° —y K~e + v e and D° —> 
TT~e + u e decays for each q 2 bin, respectively. Fits to these 


TABLE VI. Summary of the range of each q 2 bin, the num¬ 
ber of the observed events A/ 0 b S erved, the number of produced 
events N pTO du.ced , and the partial decay rate AT in each q 2 bin 
for D° —>• 7r e + ^ e decays. 


q 2 (GeV 2 /c 4 ) 

N observed 

N produced 

AT (ns" 1 ) 

(0.0, 

0.2) 

814.4 ±30.9 

1066.9 ±43.2 

0.931 ±0.038 

(0.2, 

0.4) 

697.2 ± 28.7 

935.1 ±42.8 

0.816 ±0.037 

(0.4, 

0.6) 

634.6 ± 27.7 

836.6 ±41.3 

0.730 ± 0.036 

(0.6, 

0.8) 

654.6 ± 27.8 

850.1 ±40.6 

0.742 ± 0.035 

(0.8, 

1.0) 

643.2 ± 27.3 

840.2 ± 39.9 

0.733 ± 0.035 

(1.0, 

1.2) 

578.6 ±26.3 

744.6 ± 37.7 

0.650 ± 0.033 

(1.2, 

1.4) 

509.9 ±24.7 

651.1 ±35.1 

0.568 ± 0.031 

(1.4, 

1.6) 

438.6 ± 23.2 

551.6 ±32.8 

0.481 ± 0.029 

(1.6, 

1.8) 

412.6 ±22.3 

534.7 ±31.7 

0.467 ± 0.028 

(1.8, 

2.0) 

320.9 ± 19.8 

420.6 ± 28.6 

0.367 ± 0.025 

(2.0, 

2.2) 

245.8 ± 17.0 

324.0 ± 24.7 

0.283 ± 0.022 

(2.2, 

2.4) 

165.4 ± 14.1 

229.9 ±21.7 

0.201 ±0.019 

(2.4, 

2.6) 

93.6 ± 10.7 

129.2 ± 16.7 

0.113 ±0.015 

(2.6, 

9max ) 

75.8 ± 10.0 

107.2 ±15.0 

0.094 ±0.013 
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distributions that follow the procedure described in Sec¬ 
tion IV C give the signal yields lV 0 b se rved for each q 2 bin. 
In these figures, the blue solid lines show the best fit to 
the data, while the red dashed lines show the background. 
In these fits, the background normalizations are left free. 

To account for detection efficiency and detector resolu¬ 
tion, the number of events A® bserved observed in the 4th 
q 2 bin is extracted from the relation 

-^bins 

^observed E e *4 "^produced ’ (VI-4) 

4=1 


where is the overall efficiency matrix that describes 
the efficiency and migration across q 2 bins. 

The efficiency matrix element is obtained by 


e ij 


^reconstructed ^ 


Nl 


generated 


e tag 


^trk+PID 
J corr ,ij ’ 


(VI.5) 


where V™ constructed is the number of signal Monte Carlo 
events generated in the jth q 2 bin and reconstructed in 
the 4th q 2 bin, ^generated j g num ber 0 f the sig¬ 

nal Monte Carlo events which are generated in the jth q 2 
bin, e tag is the single D° tag efficiency, and /^rrl™ is 
the efficiency correction matrix for correcting the Monte 
Carlo deviations for tracking and particle identification 
efficiencies of each element of the efficiency matrix de¬ 
scribed above. 

Table VII presents the average overall efficiency matrix 
for the D° -> K~e+ v e mode. To produce this average 
overall efficiency matrix, we combine the efficiency ma¬ 
trices for each tag mode weighted by its yields shown 
in Tab. I. The diagonal elements of the matrix give the 
overall efficiency for D° —> K~e + v e decays to be recon¬ 
structed in the correct q 2 bin in the recoil of the single D° 
tags, while the neighboring off-diagonal elements of the 
matrix give the overall efficiency for cross feed between 
different q 2 bins. Similarly, Table VIII presents the aver¬ 
age overall efficiency matrix for the D° —> TT~e + u e chan¬ 
nel. 

The number of D° —> K~(n~)e + u e semileptonic decay 
events produced with q 2 filled in the 4th q 2 bin is obtained 
from 


-^bins 

produced = E ( £ "%Ebserved, (VL6) 

j 


with a statistical error given by 


^stat 


-^bins r 

(-^produced) \ ^ ^ )r/ y^stat (-^observed )J 

\ 3 


(VI.7) 


in which er s tat(A;( bserved ) is the statistical error of 
-^observed- The partial width for the 4th bin is given by 


AT, = 


N: 


produced 


T D oN t 


tag 


(VI-8) 


where Tp >o is the lifetime of the D° meson and iV tag is the 
number of the single D° tags. 

The numbers of the signal events and g 2 -dependent 
partial widths for D° —> K~e + v e and D° —»■ n~e + v e are 
summarized in Table V and Table VI, respectively, where 
the errors are statistical only. 


B. Fitting partial decay rates to extract form 
factors 


To extract the form-factor parameters, we fit the the¬ 
oretical predictions of the rates to the measured partial 
decay rates. Taking into account the correlations of the 
measured partial decay rates among q 2 bins, the y 2 to 
be minimized is defined as 


-^bins 

^,2 ^ ^ ^^pmeasured y^ pexpected^ 1 ^ a pleasured 

i,j =1 

_ Ar expected- ) y (VL 9 ) 


where AT™ easured is the measured partial decay rate in 
4th q 2 bin, CT 1 is the inverse of the covariance matrix CV,- 
which accounts for the correlations between the measured 
partial decay rates in different q 2 bins, and Wins is the 
number of q 2 bins. The expected partial decay rate in 
the 4th q 2 bin is given by 


^expected _ 


f 1 max (») G p | Vcs(d) | , |3|j t K W f 2m2j 2 

L -24^-I/+ (« )l d Q > 

^min (4) 


(VI.IO) 

where q 2 nin r^ and q^^u) are the lower and higher bound¬ 
aries of the q 2 bin 4, respectively. In the fits, all parame¬ 
ters of the form-factor parameterizations are left free. 


We separate the covariance matrix into two parts, one 
is the statistical covariance matrix C, s ) at and the other 
is the systematic covariance matrix G ^ s . The statistical 
covariance matrix is determined by 


/^fStat _ 

— 


( i 


\Td 0 At a 


E 


^ icy. ^jo. 


(cr(V ( 


observed 


)f 


(VI.ll) 

Table IX and Table X give the statistical correlation ma¬ 
trix and relative statistical uncertainties of the measured 
partial decay rates for D° —> K~e + u e and D° —> i t~ e + v e 
decays, respectively. 

Inserting the inverse statistical covariance matrix 
(C stat )—i into Eq. (VI. 9), replacing the form factor 
f+ M (q 2 ) in Eq. (VI.IO) with different form-factor pa¬ 
rameterizations discussed in the Section II, and fitting 
to the measured partial decay rates yields the product 
°f /+ ^(0) and |I4 s ( d )| as well as the parameters of the 
form factor. 
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TABLE VII. Weighted efficiency matrix tij (in percent) for D° —> K e + v e . The column gives the true q 2 bin j , while the row 
gives the reconstructed q 2 bin i. 


6ij 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

l 

71.62 

5.39 

0.70 

0.31 

0.11 

0.03 

0.01 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

2.54 

65.14 

6.31 

0.86 

0.42 

0.15 

0.02 

0.01 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

3 

0.05 

3.19 

62.49 

6.53 

0.84 

0.37 

0.10 

0.05 

0.01 

0.02 

0.00 

0.00 

0.00 

0.01 

0.00 

0.02 

0.00 

0.00 

4 

0.02 

0.09 

3.61 

61.47 

7.07 

0.79 

0.31 

0.09 

0.01 

0.00 

0.01 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

5 

0.01 

0.02 

0.11 

3.64 

60.77 

6.93 

0.88 

0.29 

0.08 

0.01 

0.01 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

6 

0.00 

0.03 

0.05 

0.10 

3.98 

61.58 

6.81 

0.80 

0.27 

0.08 

0.01 

0.03 

0.00 

0.00 

0.01 

0.00 

0.00 

0.00 

7 

0.01 

0.01 

0.01 

0.05 

0.14 

4.16 

60.75 

6.71 

0.66 

0.27 

0.04 

0.02 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

8 

0.00 

0.00 

0.02 

0.01 

0.05 

0.16 

4.12 

60.04 

6.70 

0.75 

0.22 

0.04 

0.04 

0.02 

0.00 

0.00 

0.00 

0.00 

9 

0.00 

0.00 

0.01 

0.03 

0.01 

0.05 

0.18 

4.12 

60.61 

6.64 

0.75 

0.17 

0.03 

0.03 

0.01 

0.00 

0.00 

0.00 

10 

0.00 

0.01 

0.00 

0.00 

0.02 

0.05 

0.05 

0.17 

4.01 

59.99 

6.38 

0.62 

0.12 

0.04 

0.03 

0.00 

0.00 

0.00 

11 

0.00 

0.00 

0.01 

0.00 

0.01 

0.01 

0.04 

0.09 

0.20 

3.91 

59.97 

5.84 

0.62 

0.10 

0.01 

0.02 

0.00 

0.00 

12 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.02 

0.03 

0.05 

0.15 

3.90 

59.02 

5.61 

0.46 

0.09 

0.04 

0.00 

0.00 

13 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.03 

0.12 

0.20 

3.82 

59.59 

5.15 

0.47 

0.11 

0.03 

0.00 

14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.03 

0.10 

0.23 

3.63 

57.71 

5.06 

0.34 

0.00 

0.00 

15 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.01 

0.03 

0.07 

0.18 

3.47 

56.72 

4.57 

0.17 

0.05 

16 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.02 

0.03 

0.28 

2.82 

54.82 

3.66 

0.17 

17 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.01 

0.00 

0.02 

0.04 

0.21 

2.66 

49.34 

2.10 

18 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.01 

0.02 

0.05 

1.95 

38.30 


TABLE VIII. Weighted efficiency matrix tij (in percent) for D° —n e + v e ■ The column gives the true q 2 bin j, while the row 
gives the reconstructed q 2 bin i. 


€ij 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

l 

71.86 

4.55 

0.57 

0.04 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

2 

1.71 

67.50 

5.18 

0.47 

0.03 

0.01 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.03 

3 

0.05 

2.11 

67.78 

5.22 

0.32 

0.03 

0.00 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

4 

0.02 

0.07 

2.42 

69.05 

5.30 

0.25 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

5 

0.00 

0.03 

0.10 

2.64 

69.00 

5.11 

0.26 

0.02 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

6 

0.01 

0.02 

0.03 

0.14 

2.57 

70.20 

4.85 

0.17 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

7 

0.01 

0.02 

0.04 

0.05 

0.15 

2.67 

71.01 

4.46 

0.13 

0.01 

0.00 

0.00 

0.00 

0.00 

8 

0.01 

0.02 

0.02 

0.04 

0.10 

0.2.3 

2.75 

71.32 

4.38 

0.11 

0.00 

0.00 

0.00 

0.00 

9 

0.01 

0.01 

0.01 

0.03 

0.04 

0.12 

0.26 

2.70 

70.75 

3.73 

0.08 

0.01 

0.00 

0.00 

10 

0.00 

0.00 

0.00 

0.00 

0.03 

0.06 

0.16 

0.29 

2.55 

69.48 

3.50 

0.14 

0.02 

0.00 

11 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.07 

0.12 

0.32 

2.72 

69.08 

3.28 

0.03 

0.02 

12 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.05 

0.19 

0.38 

2.63 

65.19 

3.09 

0.01 

13 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.00 

0.01 

0.01 

0.13 

0.41 

2.49 

64.07 

2.81 

14 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.01 

0.00 

0.01 

0.02 

0.08 

0.33 

2.29 

66.78 


C. Systematic uncertainties in form factor 
measurements 

1. Systematic covariance matrix 

For each source of systematic uncertainty, an IVbins x 
IVbins covariance matrix is estimated. The total system¬ 
atic covariance matrix is obtained by summing all these 
matrices. 


(1) Number of D° tags 

The uncertainties associated with the number of the 
single D° tags are fully correlated across all q 2 bins. 
The systematic covariance contributed from the un¬ 
certainty in the number of single D° tags is calculated 


by 

^ ys (7V tag ) = ATiAT^^M) , (VI.12) 

where cr(V ta g)/A r t ag is the relative uncertainty of the 
number of the single D° tags. 

(2) D° lifetime 

The uncertainty associated with the lifetime of the 
D° meson are fully correlated across all q 2 bins, so 
the systematic covariance is calculated by 

Cij s ( T D°) = Ar * Ar j > (VLi3) 

where ct(t£ >o) is the uncertainty of the D° lifetime 
taken from PDG [6]. 
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TABLE IX. Statistical correlation matrix and relative statistical uncertainty of the measured partial decay rate in each q 2 bin 
for D° K~e + v e . 


q 2 bin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

correlation pij 

1.000 

-0.115 

1.000 










0.003 

-0.146 

1.000 









-0.003 

0.005 

-0.159 

1.000 








-0.001 

-0.005 

0.007 

-0.171 

1.000 







0.000 

-0.001 

-0.005 

0.011 

-0.172 

1.000 






0.000 

0.000 

0.000 

-0.004 

0.009 

-0.175 

1.000 





0.000 

0.000 

-0.001 

0.000 

-0.004 

0.010 

-0.174 

1.000 




0.000 

0.000 

0.000 

0.000 

0.000 

-0.004 

0.011 

-0.174 

1.000 



0.000 

0.000 

0.000 

0.000 

0.000 

-0.001 

-0.004 

0.010 

-0.171 

1.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.003 

0.008 

-0.166 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.002 

0.010 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.003 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.001 

0.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.000 

-0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

stat. uncert. 

1.31% 

1.41% 

1.49% 

1.60% 

1.61% 

1.74% 

1.77% 

1.90% 

1.97% 

2.08% 

q 1 bin 

11 

12 

13 

14 

15 

16 

17 

18 



correlation pij 

1.000 

-0.159 

1.000 










0.007 

-0.154 

1.000 









-0.002 

0.007 

-0.146 

1.000 








0.000 

-0.002 

0.006 

-0.145 

1.000 







0.000 

-0.001 

-0.001 

0.003 

-0.127 

1.000 






0.000 

0.000 

-0.001 

0.000 

0.004 

-0.119 

1.000 





0.000 

0.000 

0.000 

0.000 

-0.001 

0.006 

-0.098 

1.000 



stat. uncert. 

2.20% 

2.42% 

2.60% 

2.95% 

3.39% 

3.92% 

5.11% 

9.12% 




TABLE X. Statistical correlation matrix and relative statistical uncertainty of the measured partial decay rate in each q 2 bin 
for D° -+n~e + v e . 


q 2 bin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

correlation pij 

1.000 

-0.087 

1.000 










-0.001 

-0.105 

1.000 









0.000 

0.002 

-0.110 

1.000 








0.000 

0.000 

0.004 

-0.114 

1.000 







0.000 

0.000 

-0.001 

0.004 

-0.108 

1.000 






0.000 

0.000 

0.000 

-0.001 

0.003 

-0.104 

1.000 





0.000 

0.000 

0.000 

0.000 

-0.001 

0.003 

-0.099 

1.000 




0.000 

0.000 

0.000 

0.000 

0.000 

-0.002 

0.002 

-0.098 

1.000 



0.000 

0.000 

0.000 

0.000 

0.000 

-0.001 

-0.002 

0.001 

-0.087 

1.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.001 

-0.002 

0.000 

-0.088 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.001 

-0.004 

-0.003 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.003 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

stat. uncert. 

4.05% 

4.58% 

4.93% 

4.77% 

4.75% 

5.07% 

5.39% 

5.94% 

5.93% 

6.80% 

q 2 bin 

11 

12 

13 

14 







correlation pij 

1.000 

-0.087 

1.000 










-0.003 

-0.086 

1.000 









-0.001 

-0.002 

-0.077 

1.000 







stat. uncert. 

7.64% 

9.45% 

12.96% 

13.95% 
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(3) Monte Carlo statistics 

The systematic uncertainties and correlations in q 2 
bins due to the limited size of the Monte Carlo sam¬ 
ples used to determine the efficiency matrices are cal¬ 
culated by 

Cf S (MC Stat.) = ) £ (Observed 

\T D oN tag J ^ 

X <se r vedCov(e-J,e-;)), (VI. 14) 

where the covariance of the inverse efficiency matrix 
elements are given by [35] 

Cov^e-, 1 ) = ( VL15 ) 

ij 

(4) Form factor structure 

In order to estimate the systematic uncertainty as¬ 
sociated with the form factor used to generate sig¬ 
nal events in the Monte Carlo simulation, we re¬ 
weight the signal Monte Carlo events so that the 
q 2 spectra agree with the measured spectra. We 
then re-calculate the partial decay rates with the new 
efficiency matrices which are determined using the 
weighted Monte Carlo events. The covariance ma¬ 
trix due to this source is assigned via 

Cg"(F.F.) = <5(Ar i )<5(Ar j ), (VI. 16) 

where d(AFj) denotes the change in the measured 
partial rate in the ith q 2 bin. 

(5) Fly,max Cut 

We assign systematic uncertainties of 0.10% due to 
the Fly,max requirement on the selected events in each 
q 2 bin, and assume that they are fully correlated be¬ 
tween q 2 bins. The systematic covariance due to this 
requirement can be obtained by 

Cqf(F? 7 , m ax) = a(Ar i )a(Ar i ), (VI.17) 

where cr(Ar^) = 0.10% x ATj. 

(6) U miss fits 

The technique of fitting U m ; ss distributions affects the 
numbers of signal events observed in q 2 bins. The 
covariance matrix due to the U m i ss fits is determined 

by 

x E e i^ e Ja (*a *) 2 > (VI.18) 

Ot 

where is the systematic uncertainty of the num¬ 
ber of the signal events observed in the bin a due to 
fitting Umiss distribution, evaluated as described in 
Sect. VC5. 


(7) Tracking and PID efficiencies 

The covariance matrices for the systematic uncertain¬ 
ties associated with the tracking efficiencies and the 
particle identification efficiencies for the charged par¬ 
ticles are obtained in the following way. We first vary 
the correction coefficients for tracking (PID) efficien¬ 
cies by ±ltr, then remeasure the partial decay rates 
using the efficiency matrices obtained from the re¬ 
corrected signal Monte Carlo events. The covariance 
matrix due to this source is assigned via 

(7 8 J 8 (Tracking, PID) = <5(Ar i )<5(Ar j ), (VI.19) 

where 5(ATi) denotes the change in the measured 
partial decay rate in the ith q 2 bin. 

(8) FSR recovery 

To estimate the systematic covariance matrix associ¬ 
ated with the FSR recovery of the positron momen¬ 
tum, we remeasure the partial decay rates without 
the FSR recovery. The covariance matrix due to this 
source is assigned via 

c^ 8 (fsr) = <y(Ar i )j(Ar j ), (vi. 20 ) 

where <5 (AT,;) denotes the change in the measured 
partial decay rate in ith q 2 bin. 

(9) Single tag cancelation 

We take the systematic uncertainties associated with 
single tag cancelation as 0.12% in each q 2 bin, and 
assume they arc fully correlated between different q 2 
bins. 

The total systematic correlation matrix and relative 
systematic uncertainties for measurements of the par¬ 
tial decay rates of the two semileptonic decays of D° —> 
K~e + v e and D° —> TT~e + v e are presented in the Table XI 
and Table XII, respectively. 


2. Systematic uncertainty in measurements of form factor 
parameters 

To obtain the systematic uncertainty of the parameters 
of the form factors obtained from the fits, we add the 
matrix elements of the statistical covariance matrix and 
systematic covariance matrix together. We then repeat 
the fits to the partial decay rates. 

The central values of the form factor parameters are 
taken from the results obtained by fitting the data with 
the combined statistical and systematic covariance ma¬ 
trix together. The quadrature difference between the un¬ 
certainties of the fit parameters obtained from the fits 
with the combined covariance matrix and the uncertain¬ 
ties of the fit parameters obtained from the fits with only 
the statistical covariance matrix is taken as the system¬ 
atic error of the measured form factor parameter. 
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TABLE XI. Systematic correlation matrix and relative systematic uncertainty of the measured partial decay rate in each g 2 
bin for D° K~e+u e . 


g 2 bin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

correlation pij 

1.000 











0.284 

1.000 










0.356 

0.407 

1.000 









0.350 

0.482 

0.406 

1.000 








0.354 

0.481 

0.496 

0.412 

1.000 







0.350 

0.477 

0.485 

0.499 

0.419 

1.000 






0.353 

0.481 

0.492 

0.496 

0.513 

0.419 

1.000 





0.350 

0.477 

0.488 

0.496 

0.506 

0.509 

0.426 

1.000 




0.340 

0.465 

0.477 

0.485 

0.498 

0.493 

0.506 

0.414 

1.000 



0.332 

0.454 

0.464 

0.472 

0.482 

0.479 

0.483 

0.489 

0.388 

1.000 


0.317 

0.433 

0.442 

0.449 

0.457 

0.455 

0.461 

0.458 

0.456 

0.345 


0.297 

0.406 

0.414 

0.419 

0.428 

0.425 

0.429 

0.428 

0.419 

0.413 


0.282 

0.386 

0.395 

0.400 

0.409 

0.407 

0.411 

0.410 

0.404 

0.390 


0.264 

0.361 

0.370 

0.377 

0.386 

0.384 

0.388 

0.388 

0.381 

0.369 


0.222 

0.304 

0.312 

0.317 

0.325 

0.323 

0.326 

0.326 

0.321 

0.310 


0.211 

0.288 

0.297 

0.303 

0.312 

0.310 

0.314 

0.315 

0.310 

0.299 


0.170 

0.233 

0.241 

0.248 

0.256 

0.255 

0.259 

0.259 

0.257 

0.247 


0.120 

0.163 

0.172 

0.180 

0.188 

0.187 

0.191 

0.192 

0.192 

0.183 

syst. uncert. 

1.67% 

1.21% 

1.20% 

1.24% 

1.23% 

1.24% 

1.25% 

1.25% 

1.30% 

1.32% 

q‘ bin 

11 

12 

13 

14 

15 

16 

17 

18 



correlation pij 

1.000 











0.291 

1.000 










0.377 

0.246 

1.000 









0.349 

0.331 

0.214 

1.000 








0.295 

0.273 

0.268 

0.140 

1.000 







0.283 

0.263 

0.253 

0.243 

0.101 

1.000 






0.233 

0.216 

0.208 

0.197 

0.170 

0.063 

1.000 





0.171 

0.157 

0.152 

0.146 

0.122 

0.127 

0.019 

1.000 



syst. uncert. 

1.36% 

1.42% 

1.50% 

1.62% 

1.93% 

2.04% 

2.78% 

4.66% 




TABLE XII. Systematic correlation matrix and relative systematic uncertainty of the measured partial decay rate in each q 2 
bin for D° —> Tv~e + u e . 


g 2 bin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

correlation pij 

1.000 











0.271 

1.000 










0.286 

0.200 

1.000 









0.305 

0.299 

0.183 

1.000 








0.252 

0.246 

0.220 

0.137 

1.000 







0.276 

0.272 

0.240 

0.263 

0.134 

1.000 






0.200 

0.197 

0.174 

0.187 

0.158 

0.077 

1.000 





0.206 

0.207 

0.183 

0.197 

0.163 

0.187 

0.048 

1.000 




0.174 

0.176 

0.156 

0.168 

0.140 

0.155 

0.115 

0.031 

1.000 



0.201 

0.205 

0.181 

0.195 

0.162 

0.182 

0.129 

0.143 

0.043 

1.000 


0.169 

0.173 

0.153 

0.165 

0.138 

0.154 

0.111 

0.119 

0.101 

0.043 


0.132 

0.135 

0.120 

0.129 

0.108 

0.121 

0.088 

0.094 

0.077 

0.093 


0.144 

0.147 

0.130 

0.141 

0.117 

0.133 

0.096 

0.103 

0.089 

0.101 


0.110 

0.107 

0.093 

0.102 

0.086 

0.097 

0.069 

0.074 

0.063 

0.074 

syst. uncert. 

1.58% 

1.47% 

1.65% 

1.53% 

1.85% 

1.66% 

2.28% 

2.09% 

2.42% 

2.10% 

g 2 bin 

11 

12 

13 

14 







correlation pij 

1.000 











0.002 

1.000 










0.083 

-0.020 

1.000 









0.062 

0.047 

-0.023 

1.000 







syst. uncert. 

2.43% 

3.09% 

2.98% 

4.92% 
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D. Results of form-factor measurements 


After considering the effects of the systematic uncer¬ 
tainties on the fitted parameters, we finally obtain the 
results of these fits to the partial decay rates with each 
form-factor model. The results of these fits are summa¬ 
rized in Table XIII, where the first errors are statistical 
and the second systematic. The fits to the differential 
decay rates for D° —> K~e + v e and D° —> TT~e + i' e are 
shown in Figs. 6 and 8, respectively. 

Figures 7 and 9 show the projections of fits onto f+{q 2 ) 
for the D° —> K~e + v e and D° —> TT~e + u e decays, respec¬ 
tively. In these two figures, the dots with error bars show 
the measured values of the form factors, (q 2 ), which 
are obtained with 


u 

in which 


K(n) 


fc 2 ) = 


/AT,; 


247T 3 


Aqf G%p' 3 km \V csW \ 


„3 I f 

,, Jq 2 . PK(n)\J + 

P K(*)W = - 


K W/„2'i|2j„2 


(q 2 )\ 2 dq 2 


\f+ M (qf)\ 2 (q^ ax(i) -QLn(i))' 


(VI.21) 


(VI.22) 


where \V CS \ = 0.97343 ± 0.00015 and \V cd \ = 0.22522 ± 
0.00061 are taken from the SM constraint fit [6]. In 
the calculation of p' 3 K ^(i), f+^\q 2 ) and f+^\q 2 ) are 
computed using the two parameter series parameteriza¬ 
tion with the measured parameters. 


E. Comparison of form-factor parameters in 
different parameterizations 

For the single pole model, the fits give 

= (1.921 ± 0.010 ± 0.007) GeV/c 2 , (VI.23) 

and 


for the form factors, cx^ K is expected to be ~ 1.75 and 
a BK™ expected to be ~ 1.34 [24]. Our measured val¬ 
ues of a D ^ K and cr D_> ' 7r significantly deviate from the 
values required by the modified pole model. Table XV 
presents a comparison of our measurements of these two 
parameters with those previously measured at other ex¬ 
periments and the expected values from the Lattice QCD 
calculations. 


F. Comparison of the measured with 

LQCD predictions 

Figures 10 (a) and (b) show comparisons between our 
measured form factors and those calculated in LQCD [14] 
for D° —> K~e + v e and D° —> TT~e + u e semileptonic de¬ 
cays, respectively. From these two figures we find that, 
although our measured values of the form factors f+(q 2 ) 
and /+(q 2 ) are consistent within uncertainties with the 
LQCD predictions, our measured values of the form fac¬ 
tors significantly deviate from the most probable values 
calculated in LQCD in the regions above 0.75 GeV 2 /c 4 
and 1.5 GeV 2 /c 4 for D° —>■ K~e + v e and D° —>■ Tr~e + u e 
decays, respectively. The precision of the measured 
/+ (q 2 ) and /+(q 2 ) is much higher than that of the LQCD 
calculations. 


G. Comparison of measurements of f+{ 0) and 

/+(o) 


Using the measured f+^\o)\V cs ^)\ from the two- 
parameter series expansion fits, we obtain 


fl(0)\V cd \ 

t.) ; = 0.2001 ± 0.0026 ± 0.0016, 

f?mvo\ 


(VI.27) 


where the first error is statistical and second systematic. 
With the values of |T4s(d)l from the SM constraint fit [6], 
we find 


M^ = (1.911 ± 0.012 ± 0.004) GeV/c 2 (VI.24) 

for D° —> K~e + v e and D° —> iT~e + v e decays, respec¬ 
tively. The agreement between the extracted values of 
pole mass and the expected values (M D *+) is extremely 

poor. For comparison, Table XIV lists the values of the 
pole mass M^~^ K and measured in this analysis 

and those previously measured at other experiments. 
With the modified pole model, the fits give 

a D ~* K = 0.310 ± 0.020 ± 0.013, (VI.25) 

and 

a D ^ = 0.279 ±0.035 ±0.011 (VI.26) 

for D° —> K~e + v e and D° —> ir~e + is e decays, respec¬ 
tively. In the modified pole model (BK parameterization) 


/+( 0 ) 
/f( 0) 


0.8649 ±0.0112 ±0.0073, 


(VI.28) 


where the first error is statistical and second system¬ 
atic. This measured ratio, f+(0)/f+(0) = 0.865 ± 0.013, 
is in excellent agreement with the LCSR calculation of 
/£(0)//f(0) = 0.84±0.04 [44], but the precision is higher 
than the LCSR calculation by more than a factor of 3. 

Using the /+ r (7r)(0)|V cs ( d )| values from the two- 
parameter series expansion fits and taking the values of 
\V cs (d) | from the SM constraint fit [6] as inputs, we obtain 
the form factors 


/+(0) = 0.7368 ± 0.0026 ± 0.0036 (VI.29) 

and 

fl{ 0) = 0.6372 ± 0.0080 ± 0.0044, (VI.30) 







21 




FIG. 6. Differential decay rates for D° —> K~e + u e as func- FIG. 7. Projections on f+(q 2 ) for D° K~e + v e . 

tion of q 2 . The dots with error bars show the data and the 
lines give the best fits to the data with different form-factor 
parameterizations. 
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FIG. 8. Differential decay rates for D° — > n~e + i'e as func- FIG. 9. Projections on /+(q 2 ) for D° — » n~e + i l e- 

tion of q 2 . The dots with error bars show the data and the 
lines give the best fits to the data with different form-factor 
parameterizations. 


where the first errors are statistical and the second sys- VII. EXTRACTION OF \V ca \ AND |V C d| 

tematic. 

A. Determination of |Vj; S | and \V C d\ 


Tables XVI and XVII show the comparisons of our 
measured form factors with those measured at other ex¬ 
periments, for which different form-factor parameteriza¬ 
tions and values of |I4s(d)| have been used. Our mea¬ 
surements of these two form factors are consistent within 
errors with other measurements, but with a higher pre¬ 
cision. 


Using the values for f^(^)(0)\V cs ^) | from the two- 
parameter z-series expansion fits and in conjunction with 
/£(0) = 0.747 ± 0.011 ± 0.015 [45] and /${ 0) = 0.666 ± 
0.020 ± 0.021 [46] calculated in LQCD, we obtain 

|V cs | = 0.9601 ± 0.0033 ± 0.0047 ± 0.0239 (VII. 1) 

and 

|V cfi | = 0.2155 ± 0.0027 ± 0.0014 ± 0.0094, (VII.2) 

where the first uncertainties are statistical, the second 
ones systematic, and the third ones are due to the theo- 
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TABLE XIII. Summary of results of form factor fits to the data. 


Decay mode 

D° -> K~e+u e 

D° —> -k~ e + v e 

f+ (n) mv csW \ 

0.7209 ±0.0022 ±0.0035 
0.1475 ± 0.0014 ± 0.0005 

Single pole model 

M po ie (GeV/c 2 ) 

1.921 ±0.010 ±0.007 

1.911 ±0.012 ±0.004 


X 2 /d-o.f. 

18.8/16 

20.0/12 

Decay mode 

D° -> K~e+u e 

D° —l n~e + i / e 

f+ (n) (0)\V csW \ 
0.7163 ± 0.0024 ± 0.0034 
0.1437 ±0.0017 ±0.0008 

Modified pole model 

a 

0.309 ±0.020 ±0.013 

0.279 ±0.035 ±0.011 


X 2 /d.o.f. 

20.2/16 

12.6/12 

Decay mode 

D° -> K~e+u e 

D° —> n~e + u e 

f+ (n) (0)\V csW \ 
0.7172 ± 0.0025 ± 0.0035 
0.1435 ± 0.0018 ± 0.0009 

Two-parameter series expansion 

ri 

-2.2286 ± 0.0864 ± 0.0573 
-2.0365 ±0.0807 ±0.0257 


X 2 /d.o.f. 

19.6/16 

12.8/12 

Decay mode 

D° -> K~e + u e 

D° —t TV~e + U e 

f+ W (0) \Vcw\ 
0.7195 ± 0.0035 ± 0.0041 
0.1420 ±0.0024 ±0.0010 

Three-parameter series expansion 

ri 

-2.3338 ±0.1587 ±0.0804 
-1.8432 ± 0.2212 ± 0.0690 

ri 

3.4188 ± 3.9090 ± 2.4098 
-1.3874 ± 1.4615 ±0.4680 

X 2 /d.o.f. 

19.1/15 

11.9/11 


TABLE XIV. Comparison of measurements of the pole masses 

M° 0 Z k and 


Experiment 


M°-? e K (GeV/c 2 ) 

T (GeV/c 2 ) 

Mark III [25] 

E691 [36] 

CLEO [37] 

CLEOII [38] 

E687 (Tag) [39] 
E687 (Incl) [39] 
CLEO-II [40] 
FOCUS [41] 

Belle [29] 

BABAR [30, 31] 
CLEO-c [24] 
CLEO-c [43] 

BESIII (This work) 


1-80^0 20 ± 0.25 

2.101q;2o ±0.20 

2.101q'2o ± 0-25 

2.00 ±0.12 ±0.18 

1.97j;g-22 ± 0.07 
l-87±gi“ ± 0.07 

1.89 ±0.05+o o3 

1.93 ±0.05 ±0.03 

1.82 ±0.04 ±0.03 

1.884 ±0.012 ±0.015 

1.97 ±0.03 ±0.01 

1.97 ±0.03 ±0.01 

1.921 ±0.010 ±0.007 

1.86l° Jg ±0.05 
1.911 q3q ± 0.07 
1.97 ±0.08 ±0.04 
1.906 ± 0.029 ± 0.023 
1.95 ±0.04 ±0.02 
1.87 ±0.03 ±0.01 
1.911 ±0.012 ±0.004 


TABLE XV. Comparison of measurements of the shape parameters a D ^ >K and a D>?r 

in the modified pole model. 

Theory/Experiment 

D—*K 

a 

D—>TT 

a 

LQCD [14] 

0.50 ±0.04 ±0.07 

0.44 ± 0.04 ± 0.07 

LCSR [42] 

0 07+ 0 15 
u - u/ -0.07 

0 11+ 0 - 11 
u - 11 -0.07 

FOCUS [41] 

0.28 ±0.08 ±0.07 


Belle [29] 

0.52 ±0.08 ±0.06 

0.10 ±0.21 ±0.10 

CLEO-c(281 pb _1 ) (tagged) [24] 

0.21 ±0.05 ±0.02 

0.16 ±0.10 ±0.05 

CLEO-c(281 pb _1 )(untagged) [43] 

0.21 ±0.05 ±0.03 

0.37 ±0.08 ±0.03 

CLEO-c(818 pb" 1 ) [28] 

0.30 ±0.03 ±0.01 

0.21 ±0.07 ±0.02 

BABAR [30, 31] 

0.377 ±0.023 ±0.029 

0.268 ± 0.074 ± 0.059 

BESIII (This work) 

0.309 ±0.020 ±0.013 

0.279 ±0.035 ±0.011 


TABLE XVI. Comparison of the form factor f+(0) measured at different experiments. 


Experiment 

/+( o) 

Form-factor parameterization 

BES-II [21] 

0.78 ±0.04 ±0.03 

Single pole model 

Belle [29] 

0.695 ±0.007 ±0.022 

Modified pole model 

BABAR [30] 

0.727 ± 0.007 ± 0.005 ± 0.007 

Single pole model and modified pole model 

CLEO-c [28] 

0.739 ± 0.007 ± 0.005 ± 0.000 

Three-parameter series expansion 

BESIII (This work) 

0.7368 ± 0.0026 ± 0.0036 

Two-parameter series expansion 
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retical uncertainties in the form factor calculations. 
From the measured ratio of given in 

J-\- V U -) I Vcs | 

Eq.(VI.27) together with the LCSR calculation of 
/+(0)//f (0) = 0.84 ± 0.04 [44], we determine 


\Vcd\ 


\Vcs\ 


0.238 ± 0.004 ± 0.002 ± 0.011, 


(VII.3) 


where the first error is statistical, the second one system¬ 
atic, and the third one is from LCSR normalization. 


B. Comparison of |\4| and \V c d\ 

Table XVIII and Table XIX give comparisons of our 
measured |V4| and | 14 | with those measured at other 
experiments. Our measurements of | V cs | and | 14| are of 
higher precision than previous results from both D meson 
decays and W boson decays. 




q 2 (GeV 2 /c 4 ) 


FIG. 10. Comparisons of the measured form factors (squares 
with error bars) with the LQCD calculations [14] (solid lines 
present the central values, bands present the LQCD uncer¬ 
tainties). 


Table XX gives a comparison of our measured 
|V cd |/|V cs | with the one measured by CLEO-c [28] and 
the world average calculated with \V c d\ and |V4| given in 
PDG2014 [6]. Our measurement of the ratio is in excel¬ 
lent agreement with the world average. 


VIII. SUMMARY 

In summary, by analyzing about 2.92 fb -1 data col¬ 
lected at 3.773 GeV with the BESIII detector operated 
at the BEPCII collider, the semileptonic decays of D° —> 
K~e + v e and D° —> ir~e + i / e have been studied. From a 
total of 2793317 ± 3684 single D° tags, 70727.0 ± 278.3 
D° —> K~e + v e and 6297.1 ± 86.8 D° —> n~e + v e signal 
events are observed in the system recoiling against the 
single D° tags. These yield the decay branching frac¬ 
tions 


B(D° K~e + v e ) = (3.505 ± 0.014 ± 0.033)% 

and 

B(D° -s- TT~e + u e ) = (0.295 ± 0.004 ± 0.003)%. 

Using these samples of D° —> K~e + v e and D° —> 
7r _ e+i/ e decays, we study the form factors as a function 
of the squared four-momentum transfer q 2 for these two 
decays. By fitting the partial decays rates, we obtain the 
parameter values for several different form-factor func¬ 
tions. For the physical interpretation of the shape pa¬ 
rameters in the single pole and modified pole models, 
the values of the parameters obtained from our fits sig¬ 
nificantly deviate from those expected by these models. 
This means that the data do not support the physical 
interpretation of the shape parameter in those models. 
We choose the values of f+(0)\V cs \ and /+(0)|T4d| ob¬ 
tained with the two-parameter series expansion as our 
main result. In this case, we obtain the form factors 

/+(0) = 0.7368 ± 0.0026 ± 0.0036 

and 

/+(0) = 0.6372 ± 0.0080 ± 0.0044. 

Furthermore, using the form factors calculated in recent 
LQCD calculations [45, 46], we obtain the CKM matrix 
elements 


114,| = 0.9601 ± 0.0033 ± 0.0047 ± 0.0239 

and 


1141 = 0.2155 ± 0.0027 ± 0.0014 ± 0.0094, 

where the errors are dominated by the theoretical un¬ 
certainties in the form factor calculations. Our measure¬ 
ment of the product f?(0)\V cs \ = 0.7172±0.0025±0.0035 
(/+(0)|V cd | = 0.1435 ± 0.0018 ± 0.0009) is the most pre¬ 
cise to date and would give more precise value of |I4| 
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(|V cd |) with its precision increasing to 0.6% (1.4%) when 
the uncertainty of the value of the related form factor 
calculated in LQCD can be ignored. 

Our measurements of the branching fractions, the 
form-factor parameters and the shapes of the form fac¬ 
tor f+^\q 2 ) as a function of q 2 for D° —> K~e + v e and 
D° —y 7r e + v e decays are all the most precise to date. 
These precise measurements of f+(q 2 ), /+(g 2 ), f+{ 0), 
/+ (0) and /+ (0) /(0) are in good agreement with the 
LQCD calculations of the form factors and the LCSR 
calculations of the ratio of the form factors, but have 
higher precision than those calculated in theories based 
on QCD, and therefore will allow incisive tests of any 
future theoretical calculations. 
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TABLE XVII. Comparison of the form factor /+(0) measured at different experiments. 


Experiment 

/+( o) 

Form-factor parameterization 

BES-II [21] 

0.73 ±0.14 ±0.06 

Single pole model 

Belle [29] 

0.624 ± 0.020 ± 0.003 

Modified pole model 

CLEO-c [28] 

0.666 ± 0.019 ± 0.004 ± 0.003 

Three-parameter series expansion 

BABAR [31] 

0.610 ±0.020 ±0.005 

Three-parameter series expansion 

BESIII (This work) 

0.6372 ± 0.0080 ± 0.0044 

Two-parameter series expansion 



TABLE XVIII. Comparison of \V, 

cs | measurements. 

Experiment 

\Vcs\ 

Note 

PDG2014 [6] 

PDG2014 [6] 

PDG2014 [6] 

PDG2014 [6] 

PDG2006 [47] 

BES-II [21] 

CLEO-c [28] 

BESIII (This work) 

0.986 ±0.016 

1.008 ±0.021 

0.953 ± 0.008 ± 0.024 

0.97343 ± 0.00015 

0.94t°'^ ±0.14 

1.00 ±0.05 ±0.11 

0.985 ± 0.009 ± 0.006 ± 0.103 

0.9601 ± 0.0033 ± 0.0047 ± 0.0239 

Using D+ -)• £ + u e , D u ->• I<~£ + v t and D + ->■ I< u e + v e 
Using D+ £ + V£ 

Using D° K~£+u e and D+ ->■ K°e + v e 

Global fit in the Standard Model 

Using W boson decay 

Using D° ->■ I<~e + u e 

Using D° ->■ K~e+u e and D+ -> K°e + v„ 

Using D° ->■ K~e + v e 


TABLE XIX. Comparison of \V c d\ measurements. 

Experiment 

\Vcd\ 

Note 

PDG2014 [6] 
PDG2014 [6] 
PDG2014 [6] 
PDG2014 [6] 

CLEO-c [28] 

BESIII [48] 

BABAR [31] 

BESIII (This work) 

0.225 ±0.008 

0.220 ±0.006 ±0.010 

0.230 ±0.011 

0.22522 ±0.00061 

0.234 ± 0.007 ± 0.002 ± 0.025 
0.2210 ±0.0058 ±0.0047 

0.206 ± 0.007 ± 0.009 

0.2155 ± 0.0027 ± 0.0014 ± 0.0094 

Using D u —¥ iv~£ + i/e, D + —¥ n u e + u e and neutrino interactions 
Using D° — > n~£ + ue and D + —>- n°e + v e 

Using neutrino interactions 

Global fit in the Standard Model 

Using D° —> 'R~e + v e and D + —¥ 7r°e + ^ e 

Using D + 

Using D° —> n~e + i / e 

Using D° —>• Tv~e + u e 


TABLE XX. Comparison of |Vcd|/| V ca | measurements. 

Experiment 

\Vca\/\Vc\ 

Note 

PDG2014 [6] 

CLEO-c [28] 

BESIII (This work) 

0.228 ±0.009 

0.242 ± 0.011 ± 0.004 ± 0.012 

0.238 ± 0.004 ± 0.002 ± 0.011 

Using |V cd | = 0.225 ± 0.008 and |U CS | = 0.986 ± 0.016 
Using D —> 7re + !z e and D — > Ke + v e 

Using D° — > -K~e + u e and D° I\~e + v e 






















